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Abstract 

Using the tools and methods developed in [|J limit theorems are 
proven for the linear oscillator with random coefficients. The asymp- 
totic behaviour of the moments is studied in detail. The technique 
presented in this paper can be applied to general linear systems with 
noise and is well suited for the investigation of stochastic beam dy- 
namics in accelerators. 
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1 Linear Oscillator with Noise in the Coeffi- 
cients 



Starting point of our investigation is a nondegenerate ( u ^ ) damped 
linear oscillator under the influence of noise 

x + e(j(t) + ea)x + u 2 Q (1 + —rj(t)) x = eu £(t) (1) 

V uj J 

or written as a system of two first-order differential equations 

X = OJq z 

(2) 

z = —ujq x + e ((, — 7 z — r] x) — e 2 a z 
e is small parameter \e\ < 1 . 

The e 2 proportionality of the deterministic term in the damping part is 
connected with the fact that we will discuss the dynamics on time scales 
0(l/e 2 ) (it is the minimum time scale where the stochastic effects could es- 
sentially influence the dynamics of our oscillator). If the damping will be 
weaker it will not affect the dynamics and we can neglect it, and if it will 
be stronger it will completely change the picture of the dynamics, the typ- 
ical time scales become exponentially large 0(exp(l/e a )), a > for positive 
damping and it will require other methods (see, for example 0) that are 
beyond the scope of this paper. 

Noise has been introduced in the damping part (j(t)), as a modulation 
of the frequency ujq {i](t)) and as an external driving force £(t). 

As a model of noise we shall take stochastic processes defined by the 
following scalar products 

V (t) = b(t) ■ y(t), £(*) = h(t) ■ y(t), j(t) = d{t) ■ y(t) 

with nonrandom n-dimensional vectors b, h and d which are quasiperiodic 
in t and which can be expanded into Fourier series 

b(t) = b m exp(iv m t), L m = (b m ) 



h{t) 



(h m ) 
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d(t) = d m exp(iis m t), d- m = (d m ) 

m=—oo 

with real frequencies v m satisfying the condition 



vi + v m = 43- m + I = 0. 



In the main part of this paper the vector y(t) G R n is assumed to be a 
solution of the linear system of Ito's stochastic differential equations 

dy = Ay- dt + B dw(t) (3) 

where A and B are (n x n) and (n x r) real constant matrices respectively, 
and to(t) is an r-dimensional Brownian motion, other choices for the noise 
model will be described later on. 

As smoothness properties of the vector functions b, h and d we shall 
require the convergence of the series [] 

+oo 

) < oo, p = 0,1 (4) 
We denote 

F = [p E Z : \b p \ + \h p \ + \d p \ ^ 0} 

and introduce 



E 



K. 



+ 



^ = {(p, g) e ^ x ^ : k P + v q - fc^ol ^ 0} 

Besides the smoothness condition (HI) we also require 



mm 

k £ {0,1,2,3,4} (p,q)GF k 



inf ^ \v p + v q — k Uq\ > 5 j > 0. (5) 



The condition (|) does not exclude resonances but requires them to be iso- 
lated. This can be easily changed to some kind of Diophantine conditions 



1 For a complex vector w £ C n we use the usual spherical norm \w\ = y/w ■ w — 
yjwi ■ wf + . . . + w n ■ w* and for a (n x n) matrices with complex coefficients we shall 
use the norm \M\ = va where A is the greatest eigenvalue of the matrix M*M, which is 
compatible with the spherical norm for vectors. 
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with increasing smoothness properties (^). Note that (^j) is always satis- 
fied for periodic functions (i.e. v p = p ■ v) and for finite trigonometrical 
polynomials with arbitrary frequencies. 

In this paper we will assume that all eigenvalues of the matrix A in 
(H) have negative real parts, i.e. 

Re A fe < -5 2 s < 0, k = l,...,n (6) 

From this it follows (see, for example |J) that if the initial random vector 
i/o , independent of the r-dimensional Brownian motion w(t) — w(0) for 
< t < oo , has a normal distribution with mean value (yo) = and 
covariance matrix 



2/o • y 



J exp(rA) BB T exp{TA T ) dr = D 



then the solution of (^) y(t, j/o) is a stationary, zero-mean Gaussian process, 
with covariance function 

( exp(rA) D ; r > 

P(r) = (7) 
[ D exp(rA T ) ; r < 

Although, later on we shall not restrict the initial conditions for y 
in our noise model to be equal to the above mentioned initial conditions 
generating stationary solutions of the system (|3|) 0, all results will nevertheless 
be expressed in terms of the spectral density associated with the covariance 
function (0) ^(cu) = ^ c (u) - i ^ s (cj) where f\ 

oo 

*cM = / COs(uT)p(T)dT = - - • D 







2 For simplicity we even shall take the initial condition to be a point in n-dimensional 
Euclidean space, but if one will follow the proofs of the theorems it will be clear that 
all results of this paper will be correct if we use as initial condition an arbitrary random 
vector, independent of the r-dimensional Brownian motion w(t) — w(0) for < t < oo, 
additionally assuming that some moments of yo are finite. 

3 Note that if the matrices A and B^ 1 commute we use notation ^ for the product 
AB- 1 . 
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/UJ ± 
sm(uT)p(r)dT = A2 - • D 



For further purposes let us note that independently from real u the norm 
of the matrix admits the estimate 

|*(w)| < C (8) 

where C is some positive constant whose exact value depends on 5 2 s and 



BB 



and is unimportant for us. 



2 Special Basis in the Space of Polynomials 

Often, the influence of noise in systems such as (|l|) is studied by considering 
its influence on the unperturbed invariants of motion such as energy 



1 / o 
/ — - \ x + z~ 



or functions of the energy. For our later study of arbitrary moments we 
introduce a special time dependent (non-autonomous) basis in the space of 
polynomials. 

For all nonnegative integers m, k we define 



m / _ • _ \ k 



J m , k = exp (i (m - k) uj t) f ~ y~ ) f ~ g~ ) 

It is easy to check that the functions introduced above admit the following 
properties 

a - (I + ^0 Im,k = 

r T > — T* 

J m,fe — 1 k,m 

d. I'm, m ( 2 J 

I 12 * / r \ m +' c 

e - l-^m, k\ = Im,k ' ^m,fc = J 



Representing x and z as 

x ~\~ iz x — %z 
x = — 1 — = exp(iuj t)I 0t i + exp(-uoot)I li0 



x ~\~ iz x — %z 
z = — ^ ^— = iexp(iu t)I 0l i - iexp(-iu t)Ii,o 

and using property b we can express x m ~ k ■ z k (0 < k < m) with the help 
of the binomial theorem in the form of a linear combination of the functions 

„m—k Jt 



x m ~ k -z k = (iy 



m—k k / / \ / / 

m — k \ Ik 



J2 IK" 1 ) 9 r. ) \ n )^P(i(m-2(p + q))^ t)- Ip+ q , m -(p+ q )- 



p=0q=0 



P J \ Q 



For m 7^ k I m ,k are functions with complex values. However, we can 
also use as basis real valued functions U mt k and V rrit k which are defined by 

j j Irn, k ~\~ Ik,m jj Im,k ^k,m tt 

Um,k ^ Uk,mi Vm,k ^ 'k,m- 

Note further that the functions U mt k and V mj k can be easily expressed 
through the real valued functions U m> k and V m> k 

( x+iz \ m ( x— iz \ k I ( x+iz \ k ( x—iz \ m 
fr _ I ^ j \ 2 J ^ { 2 J { 2 J 

Um, k 



( x+iz \ 171 S x—iz \ k f x+iz \ k ( x—iz \ m 

i> _ \ 2 J { 2 J V 2 )_ { 2 J 

v m, k ^ . 

which do not depend on time t with help of the following simple formula 

1 U m ^k 

\ V m ,k 





( cos((m 


- k) u) Q t) 


— sin((m 


- k) u t) 


v v m , k ) 


^ sin((m 


- k) u t) 


cos((m 


- k) u t) 
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3 Stopped Process 



Although a suitable choice for A and B in (Q) allows one to approximate 
a wide range of spectral functions (with appropriate choice of A and B one 
can obtain for the y\ component of the vector y every spectral function 
which is the ratio of two polynomials), the solution of this equation has 
the disadvantage that it also allows with positive probability arbitrary big 
excursions during finite fixed time intervals. In order to remove this effect 
and also to apply our proof technique we have to freeze and truncate the 
process. 

Let c(e) be some positive function of e defined on the set e ^ . For 
every natural m and for every point yo G R n we introduce a random value 

r e m = r £ m {y ) = inf {t > : (t, y(t)) [0, m) x {y : \y\ < c(e)}} 

where y(t) is the solution of the system (||) which with probability one 
satisfies the initial condition y(0) = yo . So with probability one for mi < 

m 2 

< rl x < rl 2 . 

Then with probability one there exists a limit (finite or infinite) when m — > oo 
of the sequence which we will denote as 

r £ (yo) = Im^T^yo). 

In other words T e (yo) is the exit time from an open ball \y\ < c(e) for 
the solution of @ starting with probability one from initial point y$. Note 
that if the matrix BB T is nondegenerate then this exit time is finite with 
probability one. 

The joint solution of the systems (§), @ (x(t), z(t), y(t)) is a Marko- 
vian diffusion process in (n + 2) -dimensional Euclidean space. Let sf = 
min{t, t 6 } . For the noise model @ for reasons which we explained above 
we shall not study the moments of the stochastic process (x(t), z(t), y(t)) , 
but the moments of the stochastic process (x(sf), z(s%), y(sf)) (stopped 
process). We shall use the time scale O (e~ 2 ) and the difference between t 
and s e t for this time scale can be estimated with the help of the following 

Theorem A: There exist positive constants a and b so that for any 
initial point ]Jq and for any positive L 

P[r £ <^j < (exp(a|y | 2 ) + a exp(-b c 2 (e)) (9) 
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Rewriting the left hand side of the inequality (|9|) in the form 
P (t £ < 4) = P I max \t - sf\ > ) 

V E 1 ) \0<t<L/e 2 J 

we see that on the time scale considered the measure of the set where t ^ s\ 
will go to zero as e: — > if c 2 (e) — ► oo faster then b~ l \og{e~ 2 ). On the 
other hand to apply the technique of our proof we require that 

limec 3 (e) = 
so that we can not allow c(e) go to infinity too fast. 

4 Asymptotic Behaviour of Moments 

Let us introduce functions ci(m, k) of integer arguments m, k > with 
the help of 

ci(m, k) = ~T \ E {( m ~ !) ^*(^o + v p ) hp ■ (bi + id t ) - 

\Ul-Up = u> 

- (k + 1) V*(2u + u p ) (bp + idp) ■ hi - 
- k **(i/ p ) (bp + idp) - hi - k ^*(uJo + Up) hp ■ (bi - id!) } + 

+ £ (i; - id;) ■ h* - k^ T (u + up) 1; • (&; - <) } ) 

m(m — 1) x „ \ r r 

c 2 {m, k) = 2^ V (o> + fp) hp ■ hi 

c 3 (m, k) = m ( m ~ U {^*(ujo + Up) hp ■ (bi - idi) + 

vi—Vp = 3a>o 



9 



+ ^*(2lu + u p ) (b p + id p ) ■ hi} 



c^m, k) 



m(m — 1) 



J2 ^*(2uj + Up) (b p + idp) ■ (bi - id, 



c 5 (m, k) = j ]T { k ^*("p) (K + idp) ■ (bi-idi) + 



+ (k + 1) ^*(2u + u p ) (bp + idp) ■ (h - id 



[m 



1) V*(2u + Up) (b p + idp) ■ (bi + idi) } 



m 



(b; - id;) ■ (v - id?) 



m k 00 r 

c 6 (m, k) = — J2 [*(<*>o + f P ) + ^*(^o + ^) 



Zip • Zip 



p=— oo 



C7(m, fc) 



m + k 



a 



+ E 



m k 



m{u p ) + **(i/ p ) (6p + idp) ■ (bp + idp) - 



p=— oo 



- ^—^{up) (bp - id p ) ■ (bp + id p ) - ^-^/*(u p ) (b p + id p ) ■ (b p - id p ) + 



+ 



> -^m(2uo Q + u p ) + —— J -V*(2uj + Up) 



(bp + idp) ■ (bp + idp 



By using (f|) and @ it is not hard to show that c/(m, k) are correctly 
defined because the series converge absolutely for every fixed values of m 
and k . 
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Now in correspondence with an arbitrary two index array a m> k and 
nonnegative integer N 

m, fc > 0, m + k < N 

we define a vector V{a niy k] N) with (N + l)(iV + 2)/2 components with the 
help of the rule 

(m + k)(m + k + 1) 
V,(a m , fc ; N) = a m>fc , / = ^ + k + 1 

This ordering corresponds to the following ordering of the elements of the 
array a mj k (take by rows) 

a o, o 

a i,o a o,i 

^2, a i, 1 a o, 2 



OiV, OJV-1, 1 OiV-2, 2 ••• dQ,N 

Consider now the system of ordinary differential equations with constant 
coefficients 

-f- V(a m . k ; N) = K N V(a m , fc ; N) (10) 
dr 



generated with the help of the rule 



d_ 

dr 



a m ,k = c 2 {m, k) a m _ 2 ,k + e 2 *(fc, m) a m)fe _ 2 + 

Ci(m, fc) a m -i,fe + C]*(fc, m) a mjfc _i + 

c 3 (m, fc) a m _ 2 ,fe+i + c 3 *(fc, m) a m+ i jfe _ 2 + 

c 4 (m, fc) a m _ 2 ,fc+2 + c 4 * (fc, m) a m+2i fc_ 2 + 

c 5 (m, fc) a m -i,fc+i + c 5 *(fc, m) a m+lifc _! + 

c 6 (m, fc) a m _ ljfe _i + c 7 (m, fc)a m , fc 



where on the right hand side of (|TT| ) we take into account only terms with 
nonnegative indices. 
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Theorem B: Let the function c(e) satisfy the condition 

limec 3 (£:) = 

Then for arbitrary initial points Xq, zq, yo , and for arbitrary nonnegative 
integer N and positive L 



lim max 

£^0 0<t<L/e 2 



M N \e 2 s e t ) V(I m , k (s £ t ); N) - V(I m , fc (0); N)) 







where the matrix A4n(t) is the fundamental matrix solution of the system 
of linear ordinary differential equations with constant coefficients 

Remark 1: For further purposes it is important to note that the state- 
ment of the theorem B can also be written in the form 



lim max 

e^O 0<t<L/e 2 







a N {e 2 sl) -V{I m ,k{sl),N) - a N {U) ■ V(/ m>fc (0), N)] 
where cIn{ t ) is an arbitrary (N + 1)(N + 2) /2-dimensional vector satisfying 

-Kn a N (12) 



dajsr l7 , T 



dr 

Remark 2 : For physical applications one can neglect the small difference 
between t and sf (see theorem A) and we have 



V(I m ,k(t), AO) w M N (eH)V(I m ,k(0),N) 



5 Nonresonant Case 

Let us now define what we mean by nonresonant. 

Definition: We shall say that there are no resonances of order m > 
if for all integers p, q such that (p, q) G F x F 

moo® ^ u p + v q 



Definition: We shall say that there are no resonances up to order m > 
if for all integers p, q such that (p, q) G F x F 

k ujq ^ v p + v q for k = 1, . . . , m 
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In the nonresonant case only the values of c^m, k) and 07(772, k) will 
be different from zero. Introduce for them special notations 

An,k = c 7 (m, k), C m , k = c 6 (m, k) 

Note that C m , k is a symmetrical function of its arguments, i.e. C m , k = Ck, m 
and it is also a real valued function i.e. C m , k = C* m and the function ^4 mi & 
satisfies An, fc = ^4fe iJn . 

For the following let us also introduce special notations for the real and 
imaginary parts of Am, k 

t <Am, k Afc, m <o "4-m, k Afc, m 

•A-m, k 2 ' &m, k ^ 

We shall call A m , k and B mi k for reasons which will become clear later diffu- 
sion coefficient and tune shift respectively. Note that A m , m = A m ,m and 
B =0 

Theorem C: Let there he no resonances up to order 4 and let the 
function c(e) satisfy the condition 

lim£c 3 (e) = 

Then for any initial points x , z , y , for any nonnegative integers m , k 
and for any positive L 



lim max 

e^O 0<i<L/e 2 



\p=0 p=0 / 







where q = min{m, k} and the functions a™' (r) are an arbitrary 
solution of the system of linear ordinary differential equations with constant 
coefficients 

1 m, k 

aa a m, k 

, •f^m, k <*o 

ar 

j m, k 

P — — A n m ' k — C n m ' k 

•s~*-m~p, k—p Up Wra— k—p+l u p— 1 



dr 



P = 
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The proof of this theorem can be obtained from the remark to the theorem 
B with help of some straightforward calculations. 

Remark 1: We would like to note that for the study of the behaviour 
of first order moments (i.e. when m + k = 1) we actually need to avoid 
resonances in theorem B up to order 2 only. 

Remark 2: The general solution of the system of differential equations 
for the coefficients a™' k has the form 

a^\ T ) = a^-expt-A^r) 
a™> k {r) = |^ fc (0)-C m _ p+1 , fc _ p+1 ^ 

exp ( ,A m — p 5 k—p 7") 

P = 1, Q 
Choosing the initial conditions 

O) = 1, a™> k (0) =0, p = 1, ...,q 

the statement of the theorem C can be rewritten in the form 

lim max 

e^O 0<t<L/e 2 



exp{-e 2 A mtk s £ t )I mtk {s £ t ) - y m>k (0) - J]a^ fc (e 2 sf)/ m _ Pifc _ p (sf) 

In the case when m ^ k we can use the real valued functions U m> k 
and V mj k instead of the complex valued I mj k ■ Due to the symmetries 
U m ,k = U k ^ m and V mjk = —V k ^ m it is enough to consider only the case 
when m > k . So we have 

Corollary 1: Let there be no resonances up to order 4 and let the 
function c(e) satisfy the condition 

limec 3 (e) = 
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Then for any initial points x , z , y , for any nonnegative integers m, k 
satisfying m > k , and for any positive L 



lim max 

e^O 0<t<L/e 2 



/ k k \ 

M™> k (e 2 s £ t ) ■ W™> fc (sf ) - Y, M™' k (0) ■ W™' k (0) 

\p=0 p=0 I 



= 



where 



^ Um—p, k— ^ 
\ ^m— p, k— p\T) J 



M™' k (r) 



( a™> k (r) -(3™' k (T) \ 
f3^ k {r) a™> k (r) ) 



and the functions a™ ,k {r) and /3™' fc (r) are an arbitrary real solution of 
the system of linear ordinary differential equations with constant coefficients 



d 

d^ 



( <' k \ 




/ m,k \ 


= R^ k 


\p7> k ) 




K (3™> k ) 



d 

d^ 



( <*p' k \ 




( a™> k \ 






= R™' k 




n 

^m—p+1, k—p+1 


\Pp' k J 




V Pp' k J 





a 



m, k 
p-1 

m, k 



am, n 
Pp-1 



P 



,k 



R 



m, k 



_J^m—p, k—p J^m—p, k—p \ 
_gm-p, k-p _Jm-p ) k-p 



For the important particular case when we do not have an external noise 
in our system, i.e. £(t) =0 (that means that we can put h{t) = and hence 
all C my k = 0) the differential equations defining the functions a™'" 1 , a™' k 
and ft™' k admit the simple solution 



\/T*(t) ) 



exp(-A m ,kr) 



t) = exp(-A mtm T) a™ ,m (0) 

/ cos(£ m , fe r) sm(B m>k r) \ / a™> k (0) \ 

v -sin(£ m , fc r) cos{B myk r) ) { (3^ k (0) ) 
15 



a^ m {r) = 0, a^\r) = 0, /^ fc (r) = 0, p^O 
Choosing initial conditions 

a™' m (0) = 1, cC>) = 1, Po' k (0) = 

we get the following 

Corollary 2: Let = and let there he no resonances of orders 2 
and 4, and let the function c(e) satisfy the condition 

limec 3 (e) = 

Then for any initial points x , z , y , for any nonnegative integers m, k 
satisfying m > k , and for any positive L 



lim max 

£^0 0<t<L/e 2 



exp (s 2 A m , m sl) r m (st)-r m (0) 
for k = m and 



= 



lim max 

£->0 0<KL/c 2 



exp(-e 2 ^ m , fc sl)M k m {sl) 



U m ,k( S t) 

V m ,k{s e t ) 



UrnM 
Kn,fe(0) 



where 



Mi(r) 



1 cos 

V 



sin 



(A^t) 



cos (A„ ife r) J 
otherwise. 



= 



m,k 



(m — k) ujq — e 2 B„ hk 



For the important case of constant vectors b , h and d the formulae 
for A mt k, B m ,k and C m) k take the simplified form 

- m + k (m — k) 2 T . , -> (m + k) 2 T , , -> -> 



16 



+ 



m + 2mk + k 



^ c (2lu ) b-b + ^ c (2lu ) d- d+ (^ s (2lu ) - ^J(2cj )) d- b 



B m , k = -(* c (0) + *](0))d-b + 



+ 



m — k 



V a (2u )b-b + V s (2u;o) d-d - (^ c (2cu ) - ^](2uj ))d-b 



C m ,k = ^ ^ c (^o) h ■ h 

6 First and Second Order Moments 

First order moments in the nonresonant case: 

Corollary CI: Let there be no resonances up to order 2 and let the 
function c(e) satisfy the condition 

limec 3 (e) = 

Then for any initial points xq, Zq, yo and for any positive L 

( x(sf) \ ( Xq 



lim max 

£^0 0<t<L/e 2 



exp (-£^1,0 sfj M(s 



\ z 



where 



M(r) 



' cos (Jiv - e 2 B li0 )rj -sin((u;o - £ 2 #i, )t) N 
sin((^ - e 2 i3i i0 )r) cos ((u - £ 2 jBi,o)t) J 



Second order moments in nonresonant case: 

Corollary C2: Let there he no resonances up to order 4 and let the 
function c(e) satisfy the condition 

limec 3 ^) = 0. 

£^0 



Then for any initial points x , z , y G R n and for any positive L 
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lim max 

e^O 0<t<L/e 2 



r(s £ t )-r -e 2 2C ia s 
for Ai i = and 



lim max 

e^O 0<t<L/e 2 



V( S D + ^)exp(-el4 M 4) 



r + 



2Ci,i 



otherwise. 



To estimate the behaviour of the remainder of the second moments we 
shall use the functions 



#2,0 = 



x 2 — z 2 



V - xz 



Corollary C3: Let there be no resonances up to order 4 and let the 
function c(e) satisfy the condition 

lim sc 3 (e) = 

Then for any initial points x , zq, yo and for any positive L 

( U 2 , (sf) \ ( U 2 , (0) \ 



lim max 

e^O 0<t<L/e 2 



where 



exp (-e 2 A 2 ,o4) M ( s 



= 



M(r) 



cos ((2cj - e 2 B 2) ) r) - sin (J2u - e 2 B 2t ) r) 
sin ((2cj - £ 2 #2, o) t) cos ((2cj - £ 2 #2, o) r) / 
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7 Comparison with White Noise Model 



As a special case we consider white noise in this chapter i.e. 

y = Ciu(t) (13) 

where C is a real constant (n x r) matrix and w(t) is an r-dimensional 
Brownian motion. Substituting (fHf ) into (H) we have 

cfx = u;o 2; (it 

(14) 

<fc = — ujq x dt — e 2 a z dt + e C T {h — zd — xb) ■ dw(t) 

As usual for the case of multiplicative noise we shall treat the system fll4]) 
as a system of Stratonovich's stochastic differential equations. 

Introduce the matrix $ = ^CC T which plays the role of the spectral 



density for the noise model (13) and define functions Q(m, k) with the help 
of 



ci(m, k) — ^ [(m - 2k - 1) $ h p ■ b t - i (m + 2k) $ h p ■ dt} , 

u t -u p = w 



c 2 (m, fc) = ^ '- J2 ® h p- h h 

v\—v-p = 2o;o 



c 3 (m, fc) = m ( m ^ J! ®h p -(bi- id^j , 



- / , \ m(m— 1) ^ _ . r\ /r •? 

c 4 (m, fc) = i— '- $ ( 6 p + zrf pJ ' I 6 ' ~ ld - 



c 5 (m, k) = — {-(m + k)$d p -di + 

+ (k - m + 1) $ b p ■ bi + i (2k + 1) $ b p ■ dA , 
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mk 00 - 
c 6 (m, k) = — $ W 

p = —00 



m + k Amk — m(m — 1) — k(k — 1) ^ 
c 7 (m, fc) = — a + ^ ^ ^ ]T ® b p- b P + 

p = —00 



Amk + m(m + 1) + k(k + 1) ^ r r m 2 - k 2 ^ -> r 

+ E $ d p ■ d P + * ^ p ' p - 

p = — 00 p = — 00 

Consider now the system of ordinary differential equations with constant 
coefficients 

^ V(a m , fc ; iV) = t N V(a m , k ; N) (f 5) 

generated with the help of the rule 

d_ 

dr Um ' k 



c 2 


>, 


*) 




-2, fc 


+ 




77V 


Am, fc-2 


+ 


Cl 


>, 


*) 


O'm- 


-1, fc 


+ 




77V 


^m, fc— 1 


+ 


c 3 


>, 


k) 


^m- 


-2, fc+1 


+ 


c 3 *(fc, 


77V 


a m+l, fc-2 


+ 


C4 


>, 


k) 


Q"m- 


-2, fc+2 


+ 


c 4 *(fc, 


ttV 


fc-2 


+ 


c 5 


>, 


k) 


®"m- 


-1, fc+1 


+ 




77V 


a m+l, fc-1 


+ 


^6 


>, 


k) 


Q>m- 


-1, fc-1 


+ 


c 7 (m, 


fc) 


"m, fc 





(16) 



where on the right hand side of (|16|) we take into account only terms with 
nonnegative indices. 

Theorem D: For any initial points Xq, Zq, yo , for any nonnegative 
integer N and for any positive L 



lim max 

e^O 0<t<L/e 2 



M^(e 2 t)V(I m , k (t);N) - V(/ m , fe (0); N) 
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where the matrix Ai N {r) is the fundamental matrix solution of the system 
of linear ordinary differential equations with constant coefficients (fHp. 

Note that in this case and for the noise model introduced below we have 
not to distinguish between s e t and t. We also mention that if we substitute 
into the expressions of ci(m, k) the matrix $ instead of the matrix ^ 
("spectral density" of white noise) we exactly get c/(m, k) . 

8 Another Noise Model 

The technique derived in this paper can be applied to a wide class of noise 
models. As a model of noise in this section we consider the stochastic pro- 
cesses represented by the following trigonometrical polynomials^ (cosine and 
sine functions with random phases) 

<i 

r){t) = Vm^xp (i (u m t + v m ■ y)) , V-m = (Vm)* 



j(t) = 7m ex P (i {"nit + V m ■ y )) , 7-m = (7m)* 

m=—q 

with real v m and v m G R n satisfying the conditions 

\n + v m \ + \v m + vi\ = m + I = 

where the integers m, I obey m, I = —q, . . . , q. 

The vector y e R n is assumed to be a solution of the following Ito's 
system 



dy = V2Bdw(t) 



4 In order not to deal with conditions similar to (||) and (||) we consider the case of a 
finite trigonometrical sum. The extension to the case of infinite series and also the proof 
of the theorem E we leave as an exercise for the interested reader. 
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where B is a real constant (n x r) matrix and w(t) is an r-dimensional 
Brownian motion. For simplicity we assume that the (n x n) matrix B B T 
is nondegenerate and | v m \ ^ for all m = —q, . . . , q (i.e. we do not have 
deterministic harmonics in our perturbation model). 

For p = —q, . . . , q we introduce real vectors u p = B T v p e R r which 
satisfy \u p \ ^ 0, and a function Q {uj, u p ) 

fl {uj, u p ) 

and define q(m, k) as follow 





\u p 


| 2 + iu 




u p \ 


4 + UJ 2 



Ci(m, k) 



m 
1 



{(m - 1) ft (uo + v p , Up) £ P (Vi - hi) ~ 



I"p+"j+<"oI + 
\+\v p +vi\=Q 



- (A; + 1) ft {2u + Up, Hp) {rip + i j p ) & - 



- k ft (u + Up, Up) £ p (rji - kQ (v p , u p ) (r) p + i%) + 



+ { m n * ( u p> Up) (Vp - hp) £i - kSV (ujq + Up, Up) Q (rjt + hi)} 

\vi-fp+u \ + 

+ \v t -Vp\=0 ) 



c 2 (m, k) = — 11 ^ ft(^o + v P , Up) 



\vp+fl+2uiQ\ + 

+\v P +vi\=0 



c 3 {m, k) = 



m{m — 1) 



v P , Up) £ P (Vi + h) + 



+ \v P +vi\=Q 
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+ Q (2w + v p , Up) (r] p + i-j p ) ^} 



64(771, k) 



m(m — 1) 



J2 ^ (2^o + v v , Up) (rj p + i%) (r]i + i^) 



+ \Vp+Vl\=0 



c 5 (m, k) 



m 



{kQ (is p , Up) {rjp + i'jp) (r]i + iy) + 



|"p+";+2<«ol + 

\+\v p +vi\=0 



+ (fc + 1) ft {2u + Up, Up 

) (Vp + hp) (Vi + ili) ~ 



- (m - 1) ft (2cj + v v , Up) (rjp + i-f p ) (777 - i'ji)} - 



mQ * ( u p> «p) ($> - *7p) (Vi + «7«) 



\vi-v p +2io \ + 

+Wl~v p \=0 



c 6 (m, k) = — Y 



p=-q 





Up 


2 


\u p 


4 + 


(vp + u ) 2 



777 + k mk JL, 
c 7 (m, fc) = — a + — 2^ 





Tip 


2 


t7 p 


4 + "* 



|t/ p + ?:7 P | 2 + 



+ 



m 2 + k 2 X 





Up 


2 


t7 p 


4 + ^ 



+ 
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m + 2mk + k I u v 1 2 2 

4 p=- 9 |« P r + + 2 Wo) 



+ ^ „ k E l-t i'7', .. 2 [Vp-Tp + Vplp) 



p=-q 





Up 


2 


\u p 


4 + "* 



. m — k JL, u p + 2co>o , . |2 

+ 1 ph q \u P \* + (vp + 2u Jo ) 2 lVp + %lpl 

Theorem E: For any initial points Xq, Zq, yo , for any nonnegative 
integer N and for any positive L 



lim max 

£^0 0<t<L/s 2 



M N L (sH)V(I m , k (ty, N) - V(I m ,k(0); N) 







where the matrix A4n(t) is the fundamental matrix solution of the system 
of linear ordinary differential equations with constant coefficients 

-j- V(a m ,k', N) = K, N V{a m k ; N) 

generated with the help of the rule ( |T6|) in which we use c/(m, k) instead of 
&i(m, k). 



9 Proof of the Theorems 

The purpose of this section is to give a detailed proof of the theorems. 

9.1 Proof of the Theorem A 

1. From the fact that all eigenvalues of the matrix A have negative 
real parts it follows that there exists a quadratic form v( y ) satisfying the 
conditions 

C x \y\ 2 < v(y) < C 2 \y\ 2 
Ay-grad^viy) < -C 3 \y\ 2 
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Here and below Cj are some positive constants the exact values of which are 
unimportant for us. 

2. Let L be the generating differential operator of the n-dimensional 
Markovian diffusion process y i.e. 

d 1 
L = — + Ay -grad^ + - BB T grad^. • grad^ (17) 

3. Introduce the function w(y) = exp(tp v(y)) for wich 

grad a w = ip w ■ grad ^ v 



d 2 w ( d 2 v , dv dv \ 

W -7\ 7\ + W ^ ' W 



dy m dy k \dy m dy k dy m dy k 

and hence 

d 2 v , dv dv 



Lw = ^Uy. gndfv + \ £ (^ T ) m , fc ( 

V m, k = 1 \ 



+ ^ a ' o ^ 



dy m dy k dy m dy. 



From this it follows that there exist constants C4 and C5 independent of 
the value of ip such that 

Lw < i> f-C 3 |y| 2 + \c± + i;C 5 \y\ 2 )w 



Taking now ip = C3/2C5 we get 

Lw < ^(-C 3 \y\ 2 + C 4 )w (18) 

4. Define the constant x as f ne maximum of the right side in inequality 
([T8| ) with respect to the variables y 

X = max (t (_C7 3 |tf + C7 4 ) «,) (19) 



Obviously one has 



< x < — C 4 exp %p 



2 — V r C 3 
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From 



and 



it follows, that the function 



w 



w + X -J ~ 1 



(20) 



will satisfy the inequality 



Lw < 



(21) 



5. Let S( = S( A ^ . From (20) and (pi] ) immediately follows that the 
stochastic process w ( sf ) is a nonnegative supermartingale, and hence 



P k < ^ < P sup|y(sf)|>c(e) < 



vi > 



P |rapifi(^) > exp(V>Cic 2 (£:))^ < ^exp(^«(jf )) + X^j exp(-^Cic 2 (e)) 

The first two inequalities in the sequence shown above are almost obvi- 
ous, and the last one follows from the property of the stochastic process 
w (sf) to be a nonnegative supermartingale. For finishing the proof take 
a = max(x, C 2 ip) and b = C\ip . 

9.2 Proof of the Theorem B 

1. The joint solution of the systems ([|), (Q) is a Markovian diffusion 
process in the (n + 2)-dimensional Euclidean space. Let L be the gener- 
ating differential operator of this stochastic process. Separating the orders 
according to e we can represent L in the form 

L = L + e L £ + e 2 L £ 2 (22) 

where the differential operators Lq , L £ , L e 2 are defined as follows 

( 9 d\ t , d d 

L = uj \ z x — + L, L e = (f - 7 z - rj x) — , L £ 2 = -a z — 

\ ox oz J oz oz 

and L is the generating differential operator of the n-dimensional Markovian 
diffusion process y given by ([Hp. 
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2. Now we wish to show that there exist functions u e m k satisfying 

Lo ^, k = Im, k (23) 
Representing the operator L £ in the form 

L £ = U - (v + 17) exp(i o; io, 1 - (v ~ ir Y) exp(-z u t) I Xy ) 
and calculating Q 

dz = 1— exp(iuj t) I m -i,k ~ 1- exp(-iu t) I m ,k-i (24) 
and taking into account property b we have 

L £ I m> k = ~ (m £ exp(i uj t) I m _ h k - k £ exp(-i uj t) J TOi fe _x + 



[ (77 + 27) - m (r) - ij)] i m , k + 

k (r) - i'y) exp(-i2u t) I m+ i,k-i - m (77 + 27) exp(i2uj t) I m -i,k+i 
Looking for the u m> k in analogous form 

u m,k = ^ f - m a i ex P(^ Im-i, k + ka\ exp(-z ^ *) im, fe-i + 

(ma* 2 - k a 2 ) I m , k + 
m a 3 exp(i 2 u t) I m _ x k+1 - k a* 3 exp(-i 2 w t) i m +i, 



5 Starting from this point it is convenient to extend the definition of the function I Pt q 
to negative indices assuming that if p < or q < then / p , 9 = . In general, after 
this extension one has to be careful with respect to the application of the property b, but 
we have not to worry about it, because the only source of lowering indices in this paper 
is differentiation and hence if the function I p . q with negative index will appear we shall 
have automatically zero multiplyer in front of it. 
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we get the system defining the unknown ai 

L a\ + iujq ai = £ 

< La 2 = rj + ry 



(25) 



L a 3 + i 2 u a 3 = r] + z'7 



Choosing a« in the form a; = di(t) ■ y , where 

+00 

= a^ p exp(i u p t) 

and taking into account that 

dcii 



p=— 00 



(26) 



reduce the system ( pop to a system of algebraic equations for the Fourier 



we 

coefficients 



A T (u; + Up) a ljP = h p 



A T (u p ) d 2 , p = bp + id p 
A t (2cj + v P )d 3tP = bp + idp 



(27) 



where we have used the notation 

A(u) = A + iuol. 



So among the characteristic roots of A we have no purely imaginary or zero 
values the matrix A(a>) is invertiable for an arbitrary real u and hence the 
system fl2~T|) has a unique solution, which can be expressed as follows 

a — T/ \ l A*(^ + u p ) t 

a ljP = A (w + i/ p ) hp = AJ AT - - - — ^ 



a 2iP = A T (z/ P ) (6 p + idp 



A T A T + v* 



- (bp + id } 
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a 3jP = A (2lu + u p ) [bp + id p 
Using the estimate 

a- t H 



A* (2w + u p 



A T A T + (2uj + Up) 2 1 



bp + idp} 



< 



P. 



which is valid for an arbitrary real u we get for 3i tP 

1 



\%p\ < 



81 



+ 



+ 



d, 



which together with (f|) guarantees the absolute convergence and the possi- 
bility of differentiating the series (|2~6 ) term by term. 
3. Calculating L £ 2 I m> k and L e u £ m k we get 





L £ 


2 Im, k 


+ 


L £ u m 


k — 




c 2 (m, 


k) I m - 


2, k 


+ 


c|(A;, 


m) I m , k-2 


+ 


Ci(m, 


k) I m - 


1, fc 


+ 




m ) Im, k-1 


+ 


c 3 (m, 


k) Im- 


2, fc+1 


+ 




m) Im+l, k-2 


+ 


c 4 (m, 


k) Ira- 


2, fc+2 


+ 




m ) Im+2, k-2 


+ 


c 5 (m, 


ki ) Ira- 


l,fc+l 


+ 




nt) Im+l, k-1 


+ 


c 6 (m, 


ki) Im.- 


1, k-1 


+ 


c 7 (m 


k) Im, k 





(28) 
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where the functions c/(m, k) are given by the following expressions 



ci(m, k) 



k d2 — ma* 2 + (k + 1) 03 ) £ + 
A; (ai + a\)(r] + ry) - (m - 1) ai (77 - ry) exp(i cj 



c 2 (m, 


fc) = 


m 
4 


(m - 


c 3 (m, 


fc) = 


m 
4 


(m - 


c 4 (m, 


fc) = 


m 
4 


(m - 


c 5 (m, 


fc) = 


m 
4 


2a - 



c 6 (m, fc) = 
C7(m, fc) = 



m — 1) a 3 (77 + 27) exp(i 4 w £) 

fc a 2 - ma* 2 + (A; + 1) a 3 ) (77 + 27) + 
(m — 1) a 3 (r] — ry)] exp(i 2 uj t) 
f [* (ai + ^) e] 

2a + (fca 2 _ m «2 + (^ + 1)°3)( ? 7 — 27) — 
2a + (ma* 2 - k a 2 + (m + 1)03) (77 + 17) 



Note that c§(m, k) = c$(k, m) and €7(171, k) = Cj(k, m) . 

4. Introduce a(nxn) matrix K(y) = y-y T with the elements kij = yiyj. 
It is easy to check, that this matrix satisfies the equation 

L K = AK + KA T + B B T 

The usefulness of this matrix for the following is connected with the fact that 
for arbitrary complex vectors a and c 



(a-y)-(c-y) = K a • c* = Kc-a 



(29) 



5. Define the (n x n) matrix-function P u = P u (y) with the help of the 
integral 



00 

= — j exp(icjr) exp(Ar)K(y) exp(A T r) dr 



(30) 
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This integral converges because all the characteristic roots of A have negative 
real parts. Introduce the new integration variable r' = r + t , where t is 
some parameter. Then (|30D becomes 



= - I exp(io;(T / - *)) exp(A(r' - t))K(y) exp(A T (r' - *)) dr' (31) 



Differentiating ([H]) with respect to t and using that due to (f50) P^ does 
not depend on t , we obtain 

(TP 

= AT - AP W - P^A 1 - iojPu = (32) 



Calculating L P w and taking into account (|32|) we get 

L P^ = AP^ + P W A T - C7(w) = -twP w + A~ - C(w) (33) 
where we have introduced the notation 

oo 

C(u) = J exp(iuT)exp(AT)BB T exp(A T r) dr, C(0) = D 
o 

For the following let us rewrite ( j33|) in the form 

L P u + iwP w = A" - C(w) (34) 

Note that for some positive constant C\ the norms of the matrices P^ and 
C (u) can be estimated uniformly with respect to real u as follows (using |||) 

I p u I < § W\ 2 , I CM I < § (35) 

6. Now we wish to show that there exist functions gi(m, k) | satisfying 

Lgi(m, k) + q(m, k) = q(m, fc), Z = 1, . . . , 7 

and these functions have continuous first and continuous first and second 
derivatives with respect to the variables t and y respectively, and these 



3 Of course, gi(m, k) like c/(m, k) are also functions oft and y 
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functions together with the above derivative are bounded with respect to the 
variable t for fixed values of m, k, y . 

We shall show this for I = 2 and the rest can be done by analogy. 

Due to ( p9|) and the reality of the vector h (that is h = h * ) we have 

c 2 (m, k) = ai £ exp(i2 u t) = 



m(m — 1) / r , _ t> \ , , m(m — 1) / rr _ r\ / n 
— - [Ka 1 -h *J exp(z 2iu t) = — - [K a x ■ hj exp(z 2 lu t) 

Substituting in the last expression the Fourier series of the vectors d\ and 
h we transform C2(m, k) into the form 

m(m — 1) ^ / -* \ 

c 2 {m, k) = 2^ {^K ai ;P ■ hi) exp(i(v p - v x + 2u )t) = 

p,l = — oo 



m(m — 1) 



^2 Kai iP -hi+ (Ka 1>p -hijexp(i(u p -ui + 2ujo)t) 



For lj introduce the matrix 

Q{u) = -C(u) - Pui 

and denote P = — P . Taking into account ([34D we have 

L (Q(u) exp(iut)) = — K exp(iut) and LP = D — K 
Choosing now 

m(m — 1) I „ 
g 2 {m, k) = i 2^ Pa hp -hi + 

I "I— "u = 2<^o 



(Q( u p -"1 + 2 ^o) a i,p ■ hij exp(i(z/ p — Vi + 2tuo)t) 
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we obtain 



Lg 2 (m, k) + c 2 (m, k) = Da 1>p -hi 

Vi—Vp = 2u>q 

which just coincides with the expression for c 2 (m, k) if we take into account 
that 



D a\ 



D 



A*(u + v v 



h, 



Due to Q55D and (^) we have 



-V*(u + v p ) hp (36) 



max < 


1 p 1 


, max 









Ci 
5? 



and hence, as it can be easily shown, the series defining the function g 2 (m, k) 
converges absolutely with 



g 2 (m, k) | < 



m(m — 1) C\ 
4 5f 



°fj \p=-oo 



E EM ( 3? ) 



i p=— oo 



The function g 2 (m, k) is a quadratic polynomial in y and so we need 
to worry about their partial derivative with respect to t only. Expressing 



E iPpJ (Qi^i) °i,p ' hi) exp(ifJL Pt it) 



dg 2 (m, k) m(m — 1) 

di = 4 

ftp, i r o 

where // Pi j = i/ p — z/; + 2uo and using the very rough estimate for \/J, Pt i\ 
\n P ,i\ < C 2 (l + \v p \){l + \vi\), C 2 = max{l, 2|o; | } 

we get that the series defining the partial derivative with respect to the 
variable t converges absolutely with 



dg 2 (m, k) 
dt 



< 



m(m — 1) C\ 
4 ~5f 



\y\ 2 + ¥f )c 2 . 



J2 (i + Mia llP | E (i + N) 



(38) 



i p=— oo 



, p=— oo 
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Note that expressions similar to (|31)| ) hold also for Da 2 ^p and Da% )P 
Da 2 , p = -^*(u p ) (b p + idp) 



Da 3jP = -^*(2lu + Up) (bp + id } 



7. Defining u £ ^ k as 



m, k — 92 
91 
93 
94 
95 
96 



m. 



k)I 



m—2, k 



+ g£(k, m) I m ,k-2 + 



m,k)I m -i, k + 9*{k, m) I m ,k-i + 

m,k)I m _ 2 ,k+i + 9s{k, m) I m +i, k-2 + 

m,k)I m _ 2 ,k+2 + g\ ( k, m) I m +2,k-2 + 

m,k)I m -l,k+l + 9s{k, m) Im+^k-l + 

m, k) I m - Xt fc _i + g 7 {m,k)I m>k 



and acting on the function 

Im, k = Im, k + £ W^., k + ^ U m, k 

by means of the operator L we have 

LI m ,k = (Lo + £Le + £ 2 Le 2 ) I 



m, k 



Lq Im, k + £ {pQ k + I'e Im, k) + 



+ Le 2 Im,k + Le u m,k) + £ Rm,k 



2 ( t e 2 
£ 1^0 U m , k 

where for the remainder R m> & we have the expression 



2 2 

-Rm, k = Le 2 U 6 k + L £ U 6 k + SL £ 2 U £ k 



(39) 



(40) 



(41) 
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Due to our construction of the functions J m> k , u £ m ^ k and k from ( fjQ|) it 
follows that 



LI, 



m, k 



e 2 (c 2 (m, k) I rn - 2 ,k + c 2 *(k, m) I m , k -2 + 

ci(m, k) I m -x,k + c*(k, m) I m ,k-i + 

c 3 (m, k) J m _ 2 ,fc+i + c 3 *(/c, m) J m+ i jfe _ 2 + 

c 4 (m, k) I m - 2 ,k+2 + c 4 *(/c, m) J m+2 ,fc_2 + 

c 5 (m, fc) / m -i ife+ i + cg(k, m) I m+ i,k-i + 

c 6 (m, &) / m -i,fe-i + c 7 (m, fc) J m , fc ) + 



(42) 



8. So 



r \ 2 



< 1 



m-f-fc 

r \ 2 



i + 1 1 



m, k 



for p + q < m + k , then for some positive constant C 3 independent of m 

2 

and the functions u e m k , u e m k defined above and R m> k can be roughly 
estimated as follows 



u 



m, k 



< C 3 (m + k)\y\(l + \I m , k |) 



u 



m, k 



< C 3 (m + k) 2 (l + \y\ 2 )(l + \I m , k \) 



R m ,k \ < C 3 (m + k) 3 (l + \y\ 3 )(l + \ J m>fe |) 



(43) 



(44) 



(45) 



9. So I J mi m I = J mi m then for some positive constant C4 independent 
of m we have from (f42~|) 



< £ 2 C 4 m 2 (1 + J ) + e 3 \ R rt 



Using the estimate (|45|) we can rewrite (|4*6| ) in the form 



LL 



< e 2 H m (1 + /„ 



(46) 
(47) 
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and for u%, m + e u%, m we have from ( pf ) and ([I 

< Qm (1 + 4 



7/ £ -I- F 1I £ 



Here 



£ m = 2mC 3 (|y| + e2m(l + \y\ 
Let 7Y m and Q m be some positive constants. Consider the function 



(l + /„ 



/ 2 ' L ™ i 

exp — e — r 



1 - eQ n 

for which after some straightforward calculations we have 

Lv m < e 2 {(n m - n m ) + e Y~lg~ ( Gm ~ Cm )) 



(48) 



(49) 



(50) 



1 + J m>m )exp -e' 



1 - 



Choosing now 



H m = m 2 (C 4 + £8mC 3 (l + c 3 (e 
£ m = 2mC 3 (c(e) + e2m(l + c 2 « 



(51) 



and assuming that e is small enough to guarantee 1 — eQ m > we get from 
(|5l| ) that L t> m < on the set |y| < c(e) and hence 



(52) 



(V m (s £ t )} < (V m (0)) 

Using that with probability one 

1 + Im,m{s e t ) < (l+eG m )(l + I m , m {s £ t )) 
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and hence with probability one 

(l-£<? m ) (1 + Im, m «) ) exp ^-£ 2 Y~Q~ M - ^ m ( S 



(0) < (l + £S W J (1 + (0)) 

we have from fl52| ) the estimate 



(1 + I m , m (4)) < ^ (1 + / m , m (0)> exp 



— \ t — 

which for the following is rewritten in the form 



max (1 + 

- L m,m (0)> (53) 

0<t<L/e 2 



where 



Dm = 7 + £ ^T exp ( ] limP m = exp (m 2 LC A 

1-eGm \l-eQ m J e^o \ 

2 

10. Denoting « mi k = u £ m k + eu e m k and introducing the new remainder 

Rm,k = R m ,k - c 2 (m, k) u m - 2 , k ~ c 2 *(&, m) u m , k -2 

5i(m, k)u m -x >k ~ c*(k, m) u mjk -i - 

c 3 (m, k) M m _ 2 ,fc+i - c£(k, m) u m+ x tk -2 - 

c 4 (m, k) u m - 2 ,k+2 - c£(k, m) u m+2 ,k-2 - 

c 5 (m, k) u m -i t k+ i - c 5 *(/c, m) Vi.n - 

c 6 (m, fc) M m _i jfc _i - c 7 (m,k)u mtk 
which admits the estimate (as follows from (f43|)-(|44])) 



R 



m, k 



< C 5 (m + k) 3 (l + s(m + k))(l + \y\ 3 )(l + \I mik \) (54) 
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we get from 

LI m ,k = £ 2 (c 2 (m, k) I m - 2 ,k + c 2 *(/c, m) I m ,k~2 + 

ci(m, k) I m -x,k + c*(k, m) I m ,k-i + 

c 3 (m, k) I m -2,k+i + c£(k, m) I m+ i yk -2 + 

c 4 (m, k) I m -2,k+2 + c 4 * (k, m) I m+2 ,k-2 + (55) 

c 5 (m, k) / m _i,fc+i + c 5 *(fc, m) / m+ i, fc _i + 

c 6 (m, fc) J m _i )fc _i + c 7 (m, &) I m> fc ) + 



e 3 R 



m, k 



Let now A/" be as in the theorem B. Using notation V( * ; iV) we can 
rewrite (|55"D in the form of the following system 



LV(I mik ; N) = e 2 JC N V(I m , k ; N) + e 3 V(R m , k ; N) (56) 

where the matrix JCn is the same as in flTDp. 

The matrix Ai^{r) is assumed to be the fundamental matrix solution 
of (10). That means that the matrix M.~^(e 2 t) satisfies 



^-M N \e 2 t) = -e 2 M- N \eH)K. N , M N \0) = I (57) 
at 

Applying the operator L to the vector M.jl(£ 2 i) V(/ m , k (t); N) and taking 
into account (|56|), (|57|) we obtain 



L (M^(e 2 t)V(I m , k (t); N)) = M N l {e 2 t) V(R m , k (t); N) (58) 
From (^) and Dynkin's formula (see, for example ||]) it follows that 
[M- N l {e 2 sl)V{I m , k {sl)- 1 N) - V(I m , k (0); N)] 



"7 



£ 3 (/ M N \6 2 r)V(R m ,k(r); N) dr 
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or, equivalently 



M^(e 2 s £ t )V(I mik (s^;N) - V(I m>k (0);N) 



e-(V(u m>k (0);N) - M N \e 2 sl) V(u m , k (s £ t ); N)) + 



+ e" 



M^(e 2 T)V(R m , k (T); N)dr 



From (W£) we obtain 

'M^(e 2 s e t )V(I m ,k(s £ t );N) - V(/ m , fe (0); N) 



max 

0<t<L/e 2 



< 



(59) 



< e ■ max 

0<r<L 



M N \r) 



max ( 2 

0<t<L/s 2 



V(u m , k (s £ t ); N) + L V(Rm,k(s E t ); N) ) (60) 



Let us define now [m] as the smallest integer which is bigger or equal to m. 
Using fl43[), fl44|), ([54]) and simple inequalities like 

\y\ < I + \y\ 2 



1 + \Im, k\ < 2 (1 + -^j,^ 

we can obtain 

{|v(w; iV)| , |V(iL, fc ; iV)|} 



m + k < N 



max 



< 



< C 6 N 5 (1 + eN) (l + |f | 3 ) (l + J [f ]_ [f - 
Taking into account (|53| ) we have from (|6"ID 



max ( 2 

0<i<L/e 2 



V(u m , fe (s?); iV) + L V{R m , k {sl)- N) 



< 



(61) 



< C 6 N 5 (1 + eN) (l + c 3 ( £ )) (2 + L)V^ (l + 1^ [ 
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which together with 

max MzHt) < exp(C 7 N 2 ) 

0<t<L " 

and fl6~0f ) gives the final estimate we need 

max (M- N 1 {e 2 sl) V(I m , k (s e t ); N) - V(/ m , fc (0); N) 

0<t<L/e 2 \ 



< 



< 



e {l + c\e)) P (1 + 1^^(0) 



where 



P = C 6 N 5 (l + eN)(2 + L)Vr^exp(C 7 N 2 ) 



(62) 



Taking the limit e — > we have from (p2|) the proof of the theorem B 
with speed of convergence ec 3 (e). 

11. To prove the remark to the theorem B we apply the operator L to 
the function a,N{s 2 t) ■ V{I m ,k($)\ N) with the result that 



(63) 



L (a N (e 2 t) ■ V(I m , k (t); N)) = e s a N (e 2 t) ■ V(R m , k (t); N) 



The rest of the proof is just repeating all the steps from the previous point 
with the usage (|63| ) instead of (|58|). 



9.3 Sketch of the Proof of the Theorem D 

1. The solution of the system ( |i~4"| ) is a Markovian diffusion process in the 
2-dimensional Euclidean space. Let L be a generating differential operator 
of this stochastic process. Separating the orders according to e we can 
represent L in the form 



L = Lq + e L £ 2 
where the differential operators L Q and L £ 2 are defined as follows 

L — (z — -x— ) 

dt I dx dz I 



(64) 
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L e 2 = ($ d- d - aj + x $ b ■ d - $ h ■ cT) — + 



+ (§h-h - 2x§h-b- 2z§h- d + 



+ x 2 $b-b + 2xz§d-b + z 2 $ d ■ d 



dz 2 



2. Now we wish to calculate L £ 2l mk . Representing the operator L £ 2 in 
the form 



L e 2 



$K • d + $6 • d — i a + i &d ■ d exp(i luq t) Jq, i + 



+ 



$6 • d + ia — i$d ■ d 



<>xp(-/ ^„ / ) ./.!.„ ) ~ + 



+ ( <3>/i • /t — 2 &h-b + i&h ■ d exp(z cj t) Jo, 



- 2 



<M ■ 6 - i&h ■ d 



$6-6 + $6 • ci 



4i + 



+ 



$6 ■ b - $cf- d + 2i$cT- 6 exp(i 2 to t) J , 2 + 



$6-6 - $d ■ d - 2i$d-b 



exp(— i 2 ujq t) I 



2, 



dz 2 



taking into account (24), property b and the expression 



d 2 J m , k mk 



dz 2 



2 Im—X, k—X 



m(m — l) k(k — 1) , , 

exp(z 2 uj t) i m _2, fc - exp(-z 2 u t) I m ,k~2 
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we obtain that L £ 2 I m k is given by the right hand side of with q(m, k) 
as follows 



ci(m, k) 
c 2 (m, k) 
c 3 (m, k) 
Ci{m, k) 
c 5 (m, k) 



in 

2 



(m — 2k — l)®h-b — i (m + 2k) <!>h ■ d exp(i cu t) 



m(m—l) 



m(m— 1) 
2 

m(m— 1) 



$h-(b - id 
§ (b + id 



exp(z 2 t^o £) 

exp(i 3 u>o t) 



i d 



exp(i 4 cjq t) 



77) 

2 



a; - (m + £;) $ d • d + (A; — m + 1) $ b • b + 



i(2fc + l)$6-d exp(i2u t) 



c 6 (m, fc) 
c 7 (m, fc) 



m+k a _|_ 4mfc-m(m-l)-fc(fc-l) ^ g g | 



4mfc+m(m+l)+fc(fc+l) 



, $ cT- d + i 



<&d-b 



3. Now, like in the proof of theorem B, we want to show that there exist 
for fixed m and k functions gi(m, k) bounded in t and satisfying 



dgi(m, k) 
dt 



+ Ci(m, k) = C;(m, k), I 



7 



We shall show it for / = 4 and the rest can be done by analogy. 

Substituting in the expression for c^m, k) the Fourier series of the vectors 
b and d we obtain 



04(771, k) 



m(m — 1) 



[$ (bp + 4) " _ di) exp(ifi pJ t) 

p, l = — oo 



c 4 (m, k) 



m(m — 1) 



/Up, i ¥= 

where /i Pi ; = z/ p — vi + 4u; • 
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So the problem will be solved if the series 



04(771, k) 



m(m — 1) 



E 



$ [6 



dp) • (6, 



exp(i fjLpjt) 



converges and can be differentiated term by term. The absolute convergence 
and differentiability can be easily shown using (f|) and with the final 
estimates 



4 (m, k) | < 



777(777 — 1) |$| 
2 71 



+ 



dg 4 (m, k) 
dt 



< 



m(m — 1) 



l$l 



E ft 



+ 



v P= — OO 



4. The rest follows the simplified version of the proof of the theorem B. 
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